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Abstract 

A strengthened form of Schur's triangularization theorem is given 
for quaternion matrices with real spectrum (for complex matrices it 
was given by Littlewood). It is used to classify projectors {A^ = A) 
and self-annihilating operators {A"^ = 0) on a quaternion unitary space 
and examples of unitarily wild systems of operators on such a space 
are presented. Littlewood's algorithm for reducing a complex matrix 
to a canonical form under unitary similarity is extended to quaternion 
matrices whose eigenvalues have geometric multiplicity 1. 
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1 Introduction and Definitions 



We denote the set of m-hj-n matrices by Mm,n(IF), where F = (D or F = IHI, 
the skew field of real quaternions with involution 

a + bi + cj + dk = a — bi — cj — dk, a, b,c,d & IR, 

and write M„ = M„ „; A* denotes the conjugate transpose; the n-hj-n upper 
triangular Jordan block with eigenvalue A is denoted by Jri(A). 

A matrix U G M„(F) is called unitary ii U*U = I. Two matrices A and 
B G M„(F) are unitarily similar (over F) if there exists a unitary U such 
that A = U*BU; they are called unitarily equivalent if there exist unitary U 
and V such that A = UBV. 

Let A be a quaternion n-by-n matrix; A G IH is a (right) eigenvalue of 
A if there exists a nonzero v G IHl" such that Av = vX. The eigenvalues are 
defined only up to similarity: Avh = vh ■ h~^Xh for each nonzero h G IHl, so 
h~^Xh is an eigenvalue of A whenever A is. Every eigenvalue A = a+bi+cj+dk 
is similar to exactly one complex number with nonnegative imaginary part, 
namely a + ^/W+(P + (Pi p!5l Lemma 2.1]; this complex number is called a 
standard eigenvalue of A. There exists a nonsingular S G M„(1H1) such that 
S~^AS is a Jordan matrix 

J = Jn, (Ai) © • ■ ■ © JnA^k), Aj = aj + hji G (D, hj > 0, (1) 

with standard eigenvalues, determined up to permutation of Jordan blocks 
ini Chapter 3]. We will assume that Ai ^ ■ ■ ■ ^ Afc with respect to the 
following ordering in (D: 

a + bi y c + di if either a > c and b = d, or b > d. (2) 

Performing the Gram-Schmidt orthogonalization on the columns of 5* gives 
a unitary matrix U = ST, where T is an upper triangular matrix with 
positive diagonal elements; this is the QR decomposition of S. Therefore, A 
is unitarily similar to an upper triangular matrix U*AU = T~^JT having the 
same diagonal as J (Schur's theorem for quaternion matrices). For a survey 
of quaternions and matrices of quaternions, see [15j. Canonical matrices 
of sesquilinear forms and pairs of hermitian forms on a quaternion vector 
space, and selfadjoint and isometric operators on a quaternion vector space 
with indefinite scalar product, are given in jl3j . 
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This article is a result of attempts of the authors to extend to quater- 
nion matrices Littlewood's algorithm [H] for reducing a complex matrix to a 
canonical form under unitary similarity. This algorithm was discussed in |2] 
and |12]; see also and the survey Littlewood's algorithm is based 
on two statements: 

(A) Strengthened Schur Theorem. Each square complex matrix A is uni- 
tarily similar to an upper triangular matrix of the form 







Fl2 


-^13 


■ F,, - 









-^23 


■ F2S 


F = 








'^3-^ri3 












■ Fs-l^s 
















(3) 



where Ai ^ ■ ■ ■ ^ and if Aj = Aj+i then the columns of are linearly 

independent; subject to the foregoing conditions, the diagonal blocks Aj/^. 
are uniquely determined hj A. If F' is any other upper triangular matrix 
that is unitarily similar to A and satisfies the foregoing conditions, then 
F' = V*FV, where V is complex unitary and V = Vi (B ■ ■ ■ (B Vg, where each 
gM„^(€).0 

(B) Singular value decomposition. Each nonzero complex matrix A is 
unitarily equivalent to a nonnegative diagonal matrix of the form 



D = aJn^ 



l-'nt-i 



0, 



G IR, ai > ■ ■ • > at-i > 0. (4) 



If U*DV 
t/i © ■ ■ ■ S 



D. 



t-i 



where U and V are complex unitary matrices, then U 
^U\V = Ui®---® Ut-1 © V, and each Ui G M... 



Littlewood's algorithm. Let A G M„((D). Littlewood's algorithm has the 
following steps. The first step is to reduce A to the form ([3]); notice that 
the diagonal blocks and all sub-diagonal blocks of F have been completely 
reduced. We restrict the set of unitary similarities to those that preserve the 



^ This formulation is not go with the algorithm: we reduce a matrix to the form 
([3]), restrict the set of admissible transformations to those that preserve all diagonal and 
sub-diagonal blocks, then the preserving them matrices have the block-diagonal form. I 
propose the following version: "where Ai ^ • • • ^ As and if Ai = A^+i then the columns of 
are linearly independent. The diagonal blocks Xilm are uniquely determined by A. 
If V~^FV = F', where V is complex unitary and F' differs from F only in over-diagonal 
blocks, then V ^ Vi (B ■ ■ ■ (B Vg, where each Vi is rii x n^." 
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block structure of ([3]), that is, to the transformations 

F ^ V*FV, V = Vi® (5) 

The second step is to take the first nonzero superdiagonal block Fij with 
respect to lexicographically ordered indices, and reduce it to the form (jlj) by 
unitary equivalence Fij i— > V*FijVj. We make an additional partition of F 
into blocks conformal with the partition of the obtained block F^j = D, and 
restrict the set of admissible transformations ([5]) to those that preserve D 
(i.e., Vi = Ui® Ut-i ® U' and Vj = Ui ® ■ ■ ■ ® Ut-i ® V). The ith 
step of the algorithm is to take the first block that changes under admissible 
transformations and reduce it by unitary similarity or equivalence to the form 
([3]) or (jl]). We restrict the set of admissible transformations to those that 
preserve the reduced part and make additional block partitions conformal to 
the block that has just been reduced. Since we have finitely many blocks, the 
process ends on a certain matrix with the property that A is unitarily 
similar to B if and only if A°° = B°°. The matrix A°° is called the canonical 
form of A with respect to unitary similarity. 

Statement (B) holds for all quaternion matrices, that is, if A is an m-by-n 
quaternion matrix, then there exist an m-hj-m unitary matrix U, an n-hj-n 
unitary matrix V, and a diagonal matrix S with nonnegative entries such 
that A = UJ:V da Theorem 7.2]. 

In Section [2] we prove statement (A) for quaternion matrices with real 
spectrum; it need not hold for quaternion matrices with nonreal eigenvalues. 
The proof is based on modified Jordan matrices [1], which we associate with 
the Weyr characteristic [101 p. 73] of a matrix. 

In SectionOwe show that Littlewood's algorithm can be applied to quater- 
nion matrices with real spectrum that reduce to the form ([3]) with s = 2; 
Littlewood's process then consists of at most two steps. This two-step Lit- 
tlewood's process can be used to obtain the canonical forms of projectors 
[A^ = A) and self-annihilating operators {A^ = 0) on a quaternion unitary 
space. A canonical form of a complex projector was given by Dokovic [3] and 
Ikramov [5]; Ikramov's proof is based on Littlewood's algorithm. 

Unfortunately, Littlewood's algorithm cannot always be applied to 
quaternion matrices with real spectrum that reduce to the form ([3]) with 
s > 3. The reason is that in the process of reduction, one can meet a block 

^ 'Block structure' is not clear, the transformation must preserve all diagonal and 
sub-diagonal blocks. 
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with nonreal eigenvalues. The problem of classifying such matrices has the 
same complexity as the problem of classifying all quaternion matrices up to 
unitary similarity since two quaternion matrices 



Ma 




and Mb 




are unitarily similar if and only if A and B are unitarily similar. That is, if 
V is unitary and V*MaV = Mb, then V = © 1/3 © V3 by statement (A) 
for quaternion matrices with real spectrum; furthermore, Vi = V2 = and 
V{AVi = B. 

Moreover, the problem of classifying quaternion matrices up to unitary 
similarity (and even the problem of classifying, up to unitary similarity, 
quaternion matrices with Jordan Canonical Form / © 2/ © ■ ■ • © tl) has 
the same complexity as the problem of classifying an arbitrary system of 
linear mappings on quaternion unitary spaces. For example, the problem of 
classifying systems of four linear mappings 



[U, V, W are arbitrary quaternion unitary spaces) is the canonical form prob- 
lem for quaternion matrices of the form 



under unitary similarity. Indeed, by statement (A) for quaternion ma- 
trices with real spectrum, if V is unitary then V*M{A, B,C, D)V = 
M{A', B', C", D') implies V = Vi®---®V5. It also follows that Vi = V2 = 1^3, 
hence {A, B, C, D) and {A', B', C, D') are the matrices of the same system of 
linear mappings {A, B, C, V) in different orthogonal bases of U, V, W; com- 
pare with [H Sect. 2.3]. 



V 




M{A,B,C,D) 



51 I A C B 

4/ / 

3/00 

21 D 

/ 
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In particular, the problem of classifying quaternion matrices up to unitary 
similarity is equivalent to the problem of classifying m-tuples of quaternion 
matrices up to simultaneous unitary similarity 

(Ai, . . . , A J [V-^A^^, V-'A^V). 

The case for m-tuples of complex matrices was proved in [7j. Other exam- 
ples of classification problems that have the same complexity as classifying 
arbitrary systems of linear operators on unitary spaces are given in Section 

El 

In Section m we prove statement (A) for nonderogatory quaternion matri- 
ces — those matrices all of whose eigenvalues have geometric multiplicity 1 [U 
Section 1.4.4]. We then extend Littlewood's algorithm to such matrices. We 
also study the structure of their canonical matrices. 

2 A Strengthened Schur Theorem for 
Quaternion Matrices with Real Spectrum 

In this section we prove the following theorem. 

Theorem 2.1. Let A be a given square quaternion matrix and suppose that 
A has only real eigenvalues. 

(a) Then there exists a quaternion unitary matrix U such that F = U*AU 
has the form where Ai > A2 > ■ ■ ■ > are real numbers; when Aj = Aj+i, 
then rii > rij+i, and -Fj j+i is an upper triangular matrix whose diagonal 
entries are positive real numbers. 

(b) The diagonal blocks Fa = XJui are uniquely determined. The off- 
diagonal blocks Fij are determined up to the following equivalence. If V is a 
quaternion unitary matrix, then F' = V*FV has the form ^ with F/j = F^ 
(and without conditions on -F/j+ij if and only if V has the form 

V = Vi®V2®---®Vs, 

where each Vi has size rii-by-rii. 

The matrix ([3]) is a unitary variant of a modified Jordan matrix, which 
was proposed by Belitskii [1] and is obtained from the Jordan matrix by a 
simultaneous permutation of rows and columns. We define it through the 
Weyr characteristic of a matrix. 
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A list of positive integers (mi,m2, . . . ,1^^) is said to be decreasingly or- 
dered if mi > m2 > ■ ■ ■ > rrik- Given a decreasing list (mi, m2, . . . , m^), its 
conjugate is the decreasingly ordered list (ri, r2, . . . , r^) in which s = mi and 
Tj is the number of m^-'s larger than or equal to i. 

The Jordan Canonical Form Jmi(O) © ^m2(0) © ■ ■ ■ © Jmki^) of a nilpo- 
tent matrix A can be arranged so that the sizes of its Jordan blocks form 
a decreasingly ordered list m = (mi, m2, . . . ,mfc), which is called the Segre 
characteristic of A; its conjugate r = (ri, r2, . . . , r^) is called the Weyr char- 
acteristic of A [lOl p. 73]. Notice that rank(A') = ri+i + - ■ ■ + rs for 1 < / < s. 



Lemma 2.1. Let A = J,n^{0) © J„^2(0) © ■ ■ ■ © 
pose that mi > m2 > ■■• > m^. Lei (ri,r2, 
(mi, m2, . . . , mfc) . Then A is similar to 



B 



J, 



m-fe 

r 



(0) be given, and sup- 
be the conjugate of 



' Or, 


Gl2 


• 








Or, 


^23 ■ 











0,3 ■ 


















■ 









where Gi^i+i = 

Proof. Notice that 

rank(G'j j_i_iG'i+i j_i_2 



IS Ti-by^i+i. 



G, 



i+t,i+t+l) 



rank(G'. 



i+t+1- 



One checks that rank(y4') = rank(i?') for all I. It follows that A is similar to 
B. □ 



Remark 2.1. The two matrices A and B in Lemma 12.11 are permutation 
similar. To get B from A, permute the first columns of Jmi(O), Jm2(0), . . . , 
and Jm^{0) into the first k columns, then permute the corresponding rows. 
Next permute the second columns into the next columns and permute the 
corresponding rows; continue the process until B is achieved. 

Let A G M„(1H1) be given, and let J (A) be its Jordan Canonical Form 
dH). A repeated application of Lemma [2A] to the nilpotent part of J{A) — \jl 
for each of the distinct eigenvalues \j gives the following. 
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Lemma 2.2. Let A G M„(]HI) be given. Then A is similar to a unique 
matrix of the form 







Gl2 















G23 





B = 




















Gs^l,S 
















(6) 



with Ai ^ • • ■ ^ As. If Xi ^ Aj+i, then Gj^j+i = 0; otherwise, rii > nj+i and 



G 







is Hi-hy-n, 



i+i- 



Belitskii [T] called the matrix ([6]) a modified Jordan matrix and proved 
that all matrices commuting with B have an upper block-triangular form; this 
fact plays a central role in his algorithm for reducing m-tuples of complex 
matrices to a canonical form by simultaneous similarity. 

Proof of Theorem \2.1\ (a) Let A e M„(1H1) be given and suppose that 
A has only real eigenvalues, say Ai > ■ ■ ■ > A^. Lemma 12.21 guarantees 
that S~^AS = B for some nonsingular matrix 5*, and B has the form ([6]). 
Perform a Gram-Schmidt orthogonalization on the columns of S so that 
U = ST is unitary and T is an upper triangular matrix with positive diagonal 
elements. T~^ is necessarily upper triangular, and its diagonal elements are 
also positive. 



Write 








G12 

G2 





C'ls 
G23 

G3 



and 











■ 


Gs-l,S 

■■ Gs 







D12 
D2 


Di3 ■■ 
D23 ■■ 


■ Du 

■ D2s 


T = 








D3 ■ 












■■ 





Gis 
G2S 
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conformal to B. Since each Aj is real, direct computation of the product 
U*AU = T-^BT shows that 





Al/ni 


Fl2 


-^13 


■ Fu 









-^23 


■ F^s 


U*AU = T-^BT = 




















■ Fs-l^s 
















Since all the eigenvalues are real, the off-diagonal blocks satisfy 
Since T~^T = I, we have 

+ Cj^i+iDj+i = 0. 



Hence, when A,- = A 



Fi,i+i — CiGi^i+iD 



i+i- 

If Aj = Aj+i, Lemma [2.21 guarantees that > rij+i. Moreover, the form of 
Gi^i+i shows that -Fj^j+i is an upper triangular matrix whose diagonal entries 
are positive real numbers. 

(b) We now prove the uniqueness part. That the eigenvalues of F and 
their multiplicity are determined is clear. The sizes nj are also determined 
by looking at powers of F — AjJ. We can also look at a decreasingly ordered 
list of the sizes of Jordan blocks corresponding to Aj and notice that the 
conjugate of this list gives us the sizes needed. 

Let V be unitary. Suppose that F' = V*FV has the form Q and suppose 
further that F/j = F^. We claim that V is block diagonal conformal to F. 

Write V = [Vij] conformal to F (and F'). Form the products FV = VF' . 
Suppose that A^ ^ Ai. The (s, 1) block satisfies the equation 

Hence, Vgi = 0. If As_i 7^ Ai, we look at the (s — 1, 1) block and conclude 
that = 0. We proceed until Xj = Ai. 

Now, we check if A^ 7^ A2. If so, then we look at the (s, 2) block, and 
proceed as before. 
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We conclude that V is block upper triangular, but since V is unitary, V is 
block diagonal. Hence, it suffices to prove the claim when all the eigenvalues 
are the same, say A. 

As before, we write V = [Vij] conformal to F, and look at the equation 
FV = VF'. 

The (s, 1) block satisfies XVgi = AKi- However, the (s — 1, 1) block 
satisfies the equation 

XVs-1,1 + -^s-l,sKl = XVs-1,1- 

Hence, F^.i^sV^i = 0. Since -Fs_i,s is upper triangular with positive diagonal 
entries, Ki =0. 

Next, we look at the (s — 2, 1) block to get 

Vs-1,1 = XVs-2,1, 

and similarly, we conclude that = 0. Notice that the same argument 

can be used to reach the conclusion that Vi^i = for alH = 2, s. 

We then look at the (s, 2) block, (s — 1, 2) block, and so on. The conclusion 
is that V is block upper triangular. Since V is also unitary, V is in fact block 
diagonal and the sizes of the blocks in V match those of F. 



3 Applications of the Strengthened Schur 
Theorem 

A square matrix A is called a projection or idempotent ii A"^ = A; it is called 
self- annihilating a = 0. A canonical form of a complex idempotent matrix 
under unitary similarity was obtained in pl[5j, see also Section 2.3 of [14] . 

Theorem 3.1. (a) Let A be a quaternion idempotent matrix (A^ = A). 
Then A is unitarily similar to a direct sum that is uniquely determined up to 
permutation of summands of matrices of the form 

(h positive), [1] , [0] . 

(h) Let A he a self- annihilating quaternion matrix (A^ = Oj- Then A is 
unitarily similar to a direct sum that is uniquely determined up to permuta- 
tion of summands of matrices of the form 

^ (h positive), [0] . 



1 h 
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Proof, (a) Let A be a quaternion idempotent matrix. Theorem l2.1( a) ensures 
that A is similar to a matrix F that has the form ([2]). Since = A, we also 
have = F. Hence, F must have the form 



UAU* 







By Theorem 12.1( b). Fu is determined up to unitary equivalence. Now let 
-^12 = Vi'^V'z be the singular value decomposition of F12, where Vi and V2 are 
(quaternion) unitary matrices, S = diag(6i, . . . , &;)©0, and 61 >■■■>&;> 0. 
Take V = Vi(BV2 and notice that 



V*FV 



I S 




(7) 



The conclusion follows by noting that the block matrix ([7j) is permutation 
similar to a matrix that is a sum of the desired matrices. 

(b) The proof is similar to that of (a) except that A^ = means that 



F = UAU* 



and S = V1F12V2 has no zero columns. 



F12 




□ 



A self-annihilating quaternion matrix has the Jordan Canonical Form 
^2(0) © ■ ■ ■ © ^2(0) © and a simple canonical form under unitary similarity, 
as was shown in Theorem 13.1( b). What about quaternion matrices A with 
the Jordan Canonical Form J2(A) © ■ ■ ■ © J2W © A/fc? If A is real, then 
{A — Xiy = 0, and hence A is unitarily similar to a direct sum of matrices 
of the form 

„ . (6 positive), and [A] . 
U A 

However, when A ^ IR, notice that (A — A/)^ need not equal 0. Part (a) of the 
next theorem shows that the class of such matrices is unitarily wild, that is, it 
contains the problem of classifying square complex matrices up to (complex) 
unitary similarity and hence (see Section 1) it has the same complexity as 
the problem of classifying arbitrary systems of linear mappings on (complex) 
unitary spaces. Parts (b)-(d) for complex matrices were given in [SI | 
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Theorem 3.2. The problem of classifying each of the following classes of 
matrices and pairs of matrices under unitary similarity is unitarily wild: 

(a) square quaternion matrices whose Jordan Canonical Form consists 
only of Jordan blocks J2(A) and [A], where A ^ ]R is the same for all the 
matrices in the class; 

(h) square quaternion matrices A satisfying = 0; 

(c) pairs of quaternion idempotent matrices {A,B), even if A is self- 
adjoint, that is A^ = A* = A and B'^ = B; 

(d) pairs of mutually- and self- annihilating quaternion matrices {A,B), 
that is AB = BA = A^ = B^ ^ 0. 

Proof (a) Let A ^ ]R be a given eigenvalue, which we may assume is stan- 
dard, so X — X -\- yi with y > 0. To prove (a), we exhibit a mapping 
M Am e M8„(IH) such that M,N e M„(€) are (complex) unitarily 
similar if and only if Am and A^ are (quaternion) unitarily similar. 
For such a given A and M e M„(C), we define 



X 



M 



4/„ 







Mj 





3/n 


Inj 
























In 



e M4„(IH) 



and 



A/4n 




Xm 
XIau 



e M8„(1H). 



Notice that Am is similar to a direct sum of Jordan blocks J2(A). 

Now, suppose that M is unitarily similar to A^, say U*MU = N for some 
unitary [/ e M„(C). Let ^ = t/ © t/ © 17 © 17 and notice that 

{V © V)*Am{V ®V)^Am 

since jU — Uj. Hence, Am is unitarily similar to A^. 

Conversely, suppose that Am is (quaternion) unitarily similar to An, that 

is V*AmV = An for some (quaternion) unitary matrix V. We claim that 
M and N are (complex) unitarily similar. Partition the unitary matrix V 
conformal to Am, and rewrite the given condition to get 



A/4n 


Xm 




' Vu 


Vu ' 




" V^^ 


Vl2 ' 




" A/,i„ 


Xn 









. V21 


V22 . 




. V21 


V22 . 







XI An 



which yields the following equalities: 
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(i) \V^i+XmV2i^V,i\ 

(ii) \Vi2 + XmV22 = VuXn + VuX, 

(iii) XV21 = V21X, 

(iv) XV22 = V21XN + V22X. 

Writing X = x + yi, and using (iii) gives iV2i = V2ii. It follows that V21 has 
complex entries. 

Prom (i), we get y{iVii — T^i^ = —XmV2i. Write Vu — P + Qj, where P 
and Q have complex entries, so that the equality becomes 2yQk — —XmV2i. 
Write Q = [Qij] and V21 = [Aij] conformal to Xm to get 



2|/ 



Qii 


Q12 


Ql3 


Qu 


Q21 


Q22 


Q23 


Q24 


Q31 


Q32 


Q33 


Q34 




Q42 


Q43 





44 







Mj 




An 


A12 


Ai3 


Au 







Inj 






A21 


A22 


A23 










21 







A31 


A32 


A33 


^34 











In 




A41 


Ai2 


^43 


A44 



Since Aij and Qij have complex entries, we must have ^4^- = = A-^j for 
j = 1, 2, 3, 4. Equating the first two rows gives V21 = 0. Moreover, Q = as 
well, which means that Vu has only complex entries. 

Since V is unitary and since V21 = 0, we must also have V12 ~ 0. More- 
over, (iv) reduces to XV22 = V22X, so V22 has only complex entries. 

Now, (ii) reduces to XmV22 = VuXn. Write Xm = P + QmJ and 
Xn = P+QnJ, where P = diag(4/„, 3/„, 2/„, /„), and Qm and Qn have com- 
plex entries. Since Vu and V22 have complex entries, we have PV22 — VnP. 
Multiplying this equality by V*2P = PV{^ gives = (PV22){V*2P) = 
ViiP'^Vii- It follows that Vn is block diagonal; that is, it has the form 
Vn = Ci © C2 © C3 © C4. Similarly, V22 is block diagonal, and since 
PV22 = VuP, we must have V22 = Vu = Ci © C2 © C3 © C4. 

Equating the noncomplex part of (ii) gives the equality QMjV22 — 
VuQnJ- Hence, we have the following equalities: 

Cs = Ci, C3 = C2, C4 = C2, and MC4 = CiN. 
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Therefore, MCi = CiN and M is unitarily similar to A^. 

(b) Notice that using (a), the problem of classifying square quaternion 
matrices up to unitary similarity is unitarily wild. Hence it suffices to prove 
that two n X n quaternion matrices M and are unitarily similar if and 
only if the two 3n x 3n matrices 



A 



M 



In M 








-'n 





and A 



N 



M 



A% 






In 


N 








In 











0. 







are unitarily similar. One checks that A 

Suppose V* AmV = A]\r, where V is unitary. By Theorem 12.1( b). V has 
the form Vi © V2 © V3. The equality AmV = VA^ now gives Vi = V2 = V3, 
and thus V{MVi = N. 

(c) We look at the pairs of quaternion idempotent matrices 



In 





M 
M 



In-M 
In-M 



and 



In 





N In-N 

N L-N 



which are unitarily similar if and only if M and are unitarily similar, 
(d) The pairs of quaternion matrices 





' 


In ' 




' 


M ' 


^ and ^ 


' 


In ' 




' 




( 































are unitarily similar if and only if M and A^ are unitarily similar. 



□ 



4 Littlewood's Algorithm for Nonderogatory 
Matrices 

A square matrix is called nonderogatory if every (standardized) eigenvalue 
has geometric multiplicity 1, that is, its Jordan Canonical Form does not 
contain two Jordan blocks having the same standard eigenvalue [H Section 
1.4.4]. In this section, we give an algorithm for reducing a nonderogatory 
quaternion matrix A by unitary similarity to a certain matrix A°°, which has 
the property that A and B are unitarily similar if and only if A°° = B°°. 
We call such a matrix A°° the canonical form of A with respect to unitary 
similarity. 

We denote by U(JF) = {/ G F | / = /^^} the set of unitary elements of 
F, where F is M, (D, or IR. 
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Theorem 4.1. (a) Each nonderogatory quaternion matrix M is unitarily 
similar to an upper triangular matrix of the form 



A 



Ai ai2 
A2 



\i = xi + yii G €, yi> 0, 

Ai ^ ■ • ■ ^ A„ , 

ai,i+i ^ (Dj if \i = A/+1. 



(8) 







A, 



•n 



(b) The diagonal elements of A are uniquely determined. Moreover, for 
every quaternion unitary matrix S, the matrices A and A' = S*AS have the 
form ^ if and only if X[ = Ai, . . . , A'„ = A„ and 

5 = diag(si, . . . ,s„), si, . . . ,Sp eU{<C), Sp+i, . . . ,Sn el({M), (9) 

where p is such that Ai, . . . , Ap ^ IR and Xp+i, . . . , A„ G H. 

Proof, (a) The proof follows that of Theorem 12.1( a). that is, we write 
S^^MS = J, where J is the Jordan canonical form of M; and all the eigen- 
values lie in the upper half-plane. We then apply the QR factorization to S to 
obtain a unitary matrix U = ST with T an upper triangular matrix whose di- 
agonal elements are positive real numbers. The matrix A = U*MU = T~^JT 
has the desired form ([8]). Notice that when A; = A^+i, a^^z+i has the form 
—Xix + xXi + a, with a > 0. Now, if A G (D and x G M, then —Ax + xA G (Dj. 
Hence, when A; = A;+i, ai^i+i (Dj. 

(b) The proof also follows that of Theorem 12.1( b). and makes use of the 
techniques used in Theorem 13.2( a). First, notice that by the uniqueness of 
the Jordan form and the fact that the eigenvalues are ordered we must have 
X'l = Xi, . . . , X'^ = A„. 

Now, we show that S is block diagonal. If Ai 7^ A„, then the {n, 1) entries 
of AS = SA' give A„s„i = s„iAi. We express Sni = p + qj, with p, g G (D and 
conclude that Sni = since Ai and A„ are complex numbers with nonnegative 
imaginary components. Another way to look at it is that otherwise we would 
have s'^XnSni = Ai, contradicting the fact that Ai 7^ A^ and Ai and A„ are 
complex numbers with nonnegative imaginary components. 

Next, we check if Ai 7^ A„_i. If so, then we look at the (n — 1, 1) entries 
of AS = SA' to obtain A„_is„_i^i = s„_i^iAi, so that = 0. We proceed 

in this manner and conclude that Sij = whenever i > j and Aj 7^ Xj. 

It follows that S is upper block triangular, and since S is unitary, it must 
be block diagonal. Hence, it suffices to prove that the claim holds when all 
the eigenvalues coincide, that is, A = Ai = A2 = ■ ■ ■ = A^. 
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We consider two cases: A G H and A ^ H. 

Suppose A G ]R. Notice that As„i = s„iA always holds. Now, look at 
the (n — 1, 1) entries of AS = SA' to obtain As„_i,i + a„_i,„s„i = Sn_i,iA. 
Hence, a„_i^nS„i = 0, and s„i = since a„_i,„ 7^ by ([8]). Now, look at the 
{n — 2, 1) entries, then the (n — 3, 1) entries, and so on and conclude that 
Sii = whenever i > 1. 

Similarly, we look at the [n, 2) entries, the {n — 1, 2) entries, and so on 
to conclude that in fact, S is upper triangular. Since S is also unitary, S is 
also diagonal. 

Now, suppose A ^ IR. Then \ = x + yi, y > 0; the equality As„i = s„iA 
implies that isni = Snii and s„i G (D. Furthermore, As„_i.i + an-i,nSni = 
s„_i,iA implies y{isn-i,i - Sn-i,ii) + a„„i,„s„i = 0. Write s„_i,i = p + qj 
and a„_i,n = u + vj, where p,q,u,v G C to get 2yqk + uSni + vSnij = 
usni + {vsni + 2yqi)j = 0. Since usni and vsni + 'iyqi are complex numbers, 
and u (since a„_i,„ ^ (Cj by ([8])), we must have Sni = and q = (i.e.. 

Now, Asn-2,i + «n-2,n-is„-i,i = s„_2,iA implies s„_i,i = and s„_2,i G (C. 
We repeat this process until we obtain Sij = for all i > j and Sjj G (D. Since 
5" is unitary, S is diagonal. □ 



An algorithm for reducing a matrix A of the form ([8]) to canonical 
form with respect to unitary similarity 

By Theorem 14.1( b). the diagonal entries of A are uniquely determined. 
Furthermore, all unitary similarity transformations that preserve the trian- 
gular form of A and its diagonal entries have the form: 

A ^ S*AS, S ego= X ■ ■ ■ X W((D) x W(lHl) x ■ ■ • x U{m) . (10) 

p n—p 

We successively reduce the off-diagonal entries aij (z < j) to a canonical 
form in the following order: 

(11) 

On each step, we use only those transformations (|T0|) that preserve the al- 
ready reduced entries. 
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Suppose that all entries that precede air in the sequence ( ITTi) have been 
reduced, and let all the transformations (llUp that preserve the entries pre- 
ceding air have the form 

A^S*AS, S eg = {S = diag{si,...,sn) eGoin}, (12) 

where 7^ is a set of relations of the form 

Si G (D, Si G IR, Si = Sj, or Si = sj^ E (D. (13) 

We reduce air to canonical form a^^ by transformations f|T2|) and show that 
all transformations f|T2|) that preserve a^^ have the form A t— > S'*yl5', 5* e 
^' = G ^ I A7^}, where ATZ consists of relations of the form f|T3l) : this is 
required for the correctness of the induction step. 

As follows from the form of relations fll3p . for every i G {1, . . . , n} there 
exists Fj G {M, (D, IR} such that 

{si\S eg} = Ui¥i). 

If a/r is not changed by transformations (fT2|) . we set a^^ = a^^ and ATZ = 0. 
Denote by P the set of positive real numbers and suppose that 

air = Zi+ Z2j, zi,Z2e (D, (14) 

was changed by transformations (fT2l) . We have the following cases. 

1) = IHl or Fr = HI. If TZ does not imply si = Sr, then we reduce air 
to the form aj^, = s^^airSr G P and obtain ATZ = {si = s^}. If si = Sr follows 
from 71, then F; = F,. = M, take a'l^ = s~[^airSi G € with a nonnegative 
imaginary component (note that a^^ ^ P, otherwise air is not changed by 
admissible transformations) and obtain ATZ = {si G (D}. 

2) F; = Fr = (D and TZ does not imply Si = Sr or Si = s~^. Then by 

dHj) a'l^ = S^^airSr = Sl'^ZiSr + {s'[^Z2Sr)j = ZiS^^ Sr + Z2Sj^^ S'^ j. If Z1Z2 ^ 0, 

we make 2;^ G P, then si = Sr (to preserve z[), next make Z2 E F, then 
si = Sr = ±1; we obtain a'^^ = PI + Pj and ATZ = {si = G P}. 

If 7^ = ^2, "we make a'l^ G P and obtain ATZ = {si = Sr}. If zi = 0, 
then 2:2 7^ (otherwise air is not changed by admissible transformations) we 
make aj^, G Pj and obtain ATZ = {si = s~^}. 

3) F; = Fr = (D, = Sr or si = s^^ . If si = Sr, then a'^^ = zi + Z2s'['^j, 
make aj^ G (D + Pj and obtain A7^ = {s; G P}. 
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If S; = s^^, then a^^ = 2:1 ^ + Z2j, make a^^, G P + €j and obtain 
A7^ = {si G ]R}. 

4) Either = C anc? F^, = ]R, or F^ = ]R and = C. Make = 
1 + 2;2j or a'l^ = j and obtain ATI = {si = Sj. & IR}. 

5) Fi = Fr = H. Make = z[ + z'^j with z[ y (see ([2])), or a^^ = 
with Z2 >- 0, and obtain ATI = {si = Sr}. 

The process ends with the reduction of ai„. We denote the matrix ob- 
tained by it is the canonical form of A with respect to unitary similarity. 
At each step we reduced an entry to a form that is uniquely determined by 
the already reduced entries and the class of (quaternion) unitarily similar to 
A, and so we obtain the following theorem: 

Theorem 4.2. Two nonderogatory quaternion matrices A and B are uni- 
tarily similar if and only if A°^ = 5°° . 

For a canonical nxn matrix A = A°°, its graph T{A) has vertices 1, . . . ,n, 
and / and r are jointed by an edge if and only if the relations si = Sr and 
si = do not follow from the condition of preserving the entries of A that 
precede air, but one of them follows from the condition of preserving air (i-e., 
si = Sr or Si = s~^ is contained in ATZ; see the cases 1-5). Notice that there 
is an edge {i,i + 1) if Aj = Aj+i in since then aj^j+i ^ (Dj. 

A square matrix A is called unitarily indecomposable if it is not unitarily 
similar to a direct sum of square matrices. 

Theorem 4.3. (a) The graph of each canonical matrix is a union of trees. 
Any union of trees with numbered vertices can be the graph of a canonical 
matrix. 

(b) A canonical matrix is unitarily indecomposable if and only if its graph 
is a tree. Moreover, let the graph T{A) of a canonical matrix A be the union 
of m trees Fj (1 < i < m) with the vertices vn < Vi2 < ■ ■ ■ < Vi^. Rearrange 
the columns of A such that their old numbers form the sequence 

Vll , . . . , Vir-^ ; V21, ■ ■ ■ , V2r2 ! • • • ! ^ml , • • • , '^mr™, , 

then rearrange its rows in the same manner. The matrix obtained has the 
form Ai Q) ■ ■ ■ ® Am, where each Ai is a unitarily indecomposable canonical 
Ti X Tj matrix. 

Proof, (a) Let the graph r{A) of a canonical matrix A have a cycle f 1 — V2 — 
Vp—vi (p > 2), and let, say, a^^^^^}, . . . , a{^^_^^^} precede aj^^^j in the 
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sequence ( ITTi) . where a^ij} denotes aij if i < j and aji if i > j. Then the 
equahty s^^ = sf^ follows from the condition of preserving the entries of A 
that precede a^y^^^^}, a contradiction to the existence of the edge Vp — vi. 

Let a graph F with vertices 1, . . . , be a joint of trees. Take A of the form 
([8]), in which Ai = ni, X2 = {n — . . . , A„, = z, and, for every I < r, air = 1 
if there is the edge I — r and air = otherwise. Clearly, A is a canonical 
matrix and T{A) = T. 

(b) Let A be a canonical matrix. Since a^ = whenever i and j are not 
connected in T{A), the graph of a unitarily indecomposable canonical matrix 
is a tree. It follows from the algorithm of reduction to canonical form that 
if T{A) is not a connected graph, then A can be reduced to a direct sum of 
unitarily indecomposable canonical matrices by simultaneous permutation of 
its rows and columns. □ 
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